Conservative corrections to the innermost stable circular orbit (ISCO) of a Kerr black 
hole: a new gauge-invariant post-Newtonian ISCO condition, and the ISCO shift due 
to test-particle spin and the gravitational self-force 

Marc FavatfQ 

Jet Propulsion Laboratory, 4800 Oak Grove Drive, Pasadena, California 91109, US$[ and 
Theoretical Astrophysics, 350-17, California Institute of Technology, Pasadena, California 91125, USA 

(Dated: 12 October 2010) 

The innermost stable circular orbit (ISCO) delimits the transition from circular orbits to those 
that plunge into a black hole. In the test-mass limit, well-defined ISCO conditions exist for the 
Kerr and Schwarzschild spacetimes. In the finite-mass case, there are a large variety of ways to 
define an ISCO in a post-Newtonian (PN) context. Here I generalize the gauge-invariant ISCO 
condition of Blanchet and Iyer [Classical Quantum Gravity 20, 755 (2003)] to the case of spinning 
(nonprecessing) binaries. The Blanchet-Iyer ISCO condition has two desirable and unexpected 
properties: (1) it exactly reproduces the Schwarzschild ISCO in the test-mass limit, and (2) it 
accurately approximates the recently calculated shift in the Schwarzschild ISCO frequency due to 
the conservative-piece of the gravitational self-force [Barack and Sago, Phys. Rev. Lett. 102, 191101 
(2009)]. The generalization of this ISCO condition to spinning binaries has the property that it also 
exactly reproduces the Kerr ISCO in the test-mass limit (up to the order at which PN spin corrections 
are currently known). The shift in the ISCO due to the spin of the test-particle is also calculated. 
Remarkably, the gauge-invariant PN ISCO condition exactly reproduces the ISCO shift predicted by 
the Papapetrou equations for a fully-relativistic spinning particle. It is surprising that an analysis of 
the stability of the standard PN equations of motion is able (without any form of "resummation" ) 
to accurately describe strong-field effects of the Kerr spacetime. The ISCO frequency shift due to 
the conservative self-force in Kerr is also calculated from this new ISCO condition, as well as from 
the effective-one-body Hamiltonian of Barausse and Buonanno [Phys. Rev. D 81, 084024 (2010)]. 
These results serve as a useful point of comparison for future gravitational self-force calculations in 
the Kerr spacetime. 

PACS numbers: 04.25.Nx, 04.25.-g, 04.25.D-, 04.30.-w 



I. INTRODUCTION, MOTIVATION, AND 
SUMMARY 

The innermost stable circular orbit (ISCO) is a point 
of dynamical instability in black hole (BH) spacetimes 
that separates stable, circular, and bound geodesic or- 
bits from those that "plunge" into the BH event hori- 
zon. The location of the ISCO can be quantified in a 
gauge-invariant manner by specifying its orbital angular 
frequency as measured by a distant observer. For a test- 
particle in the Schwarzschild spacetime, this frequency 
occurs at m2^l = 6~ 3 / 2 , where rri2 is the mass of the BHQ 
The location of the ISCO is important in the context 
of quasicircular, inspiralling compact binaries (an im- 
portant source for ground and space-based gravitational- 
wave detectors) because it represents the point where the 
character of the orbit (and hence the gravitational waves) 
abruptly changes. Because of this, the ISCO frequency 
is often taken as the termination point of inspiral tem- 
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plates. The ISCO is also important because its location 
encodes (potentially observable) information about the 
strong-gravity region of the BH spacetime. 

What happens if we no longer have a geodesic or- 
bit? When dissipation (i.e., radiation-reaction) is in- 
cluded, the location of the ISCO is no longer precisely 
quantifiable — it becomes "blurred" into a transition re- 
gion (in orbital radius or frequency) separating the adia- 
batic inspiral from the plunge [H, 0- However, if we con- 
sider only conservative corrections to geodesic motion, a 
precise ISCO can (in some cases) continue to exist. In 
particular here we will consider two types of conserva- 
tive corrections to geodesic motion: (i) the gravitational 
self-force (GSF; a force arising from the point-particle's 
finite mass which causes it to deviate from geodesic mo- 
tion) and (ii) the force due to the spin of the test-bodyH 

Calculations of the GSF are motivated by the need 
to model extreme-mass-ratio inspirals (EMRIs), an im- 
portant source for the planned Laser Interferometer 
Space Antenna (LISA) [3] consisting of a compact ob- 
ject (mi ~ 1-IOOMq) inspiralling into a massive BH 
(m 2 ~ 10 4 -10 7 M Q ) with mass ratios q < 10 4 . Comput- 
ing the GSF is challenging (see 0-0] for reviews and refer- 
ences) , but several groups have had recent success 
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In particular, one of the concrete results to emerge from 
the self-force program has been the calculation by Barack 
and Sago (BS) of the shift in the ISCO frequency due 
to the conservative GSF in the Schwarzschild spacetime 
[Til Il5j ]. This result is especially interesting because it 
supplies a gauge-invariant, exact strong-field result that 
is only computable using the full self-force formalism. 
(This is in contrast to standard BH perturbation the- 
ory calculations, which only provide access to the time- 
averaged dissipative pieces of the self- force.) The result- 
ing conservative GSF ISCO frequency shift can be ex- 
pressed in the form 



MQ = 6~ 3/2 [l + 7?c GSF (0) + 0(r] 2 )} 



(1.1) 



where BS calculated the value c (0) = 
1.2512(±0.0004). This value can be used to com- 
pare different GSF codes, and to set constraints on 
the effective-one- body (EOB) @, [TM3 formalism (see 

In Ref. ;21i] I compared the above GSF ISCO shift 
with ~ 15 distinct post-Newtonian (PN) or EOB meth- 
ods for computing the ISCO. Among those methods, two 
approaches — based on the EOB formalism and the stan- 
dard PN equations of motion — have especially desirable 
features. In particular, the best agreement (~ 10% error) 
with the BS result was found using a version of the EOB 
formalism in which a pseudo-4PN term is added to the 
effective metric and calibrated with the Caltech/Cornell 



numerical relativity simulations [22[. This method also 
adequately predicted (with ~ 16% error) the ISCO fre- 
quency for equal-mass binaries as computed from se- 
quences of quasicircular initial data [23| . However, in the 
absence of calibration, the method which most accurately 
reproduced the BS result was the gauge- invariant ISCO 
condition of Blanchet and Iyer [24[ o This condition is de- 
rived from a stability analysis of the 3PN (nonspinning) 
equations of motion; it takes the form 



Co = 1 - 6x + Urjx 2 
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where x = (Mil) 2 / 3 , and Cq > is required for stable, 
circular orbits to exist. The ISCO is found by solving 
Co = for x (or f2). The resulting value for the conser- 
vative GSF ISCO shift was found to be (H 
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(1.3) 



3 In addition to the above-mentioned reasons, these two methods 
for computing the ISCO are also preferred over the other ap- 
proaches examined in [2lj | because: (i) the error in the ISCO 
computed via these methods decreases monotonically as the PN 
order is increased, and (ii) they each are derived from equations 
of motion that allow for a complete description of the two-body 
dynamics. 



which differs from the exact BS result by 14.7%. 

The above PN ISCO condition is especially interesting 
because it exactly reproduces the Schwarzschild ISCO 
(x = 1/6 or 7712^ = 6~ 3 / 2 ) in the test-particle limit. 
It is surprising that a condition derived from the PN 
equations of motion can reproduce a strong-field result 
like the ISCO@ For example, a standard way to compute 
the ISCO in a PN context is by finding the minimum of 
the circular-orbit energy^ In the test-mass limit, the PN 
expansion of the circular-orbit energy, 
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converges slowly: to get within 8% of the exact result 
(x = 1/6) one needs to truncate the above expression at 
4PN order or higher. 

Part of the motivation for developing "resummation" 
methods was to cure this problem while also providing a 
means to compute the ISCO for finite mass-ratio bina- 
ries. For example, Kidder, Will, and Wiseman (2(| [27j 
modified the PN equations of motion by replacing the 
0(rp) terms with the corresponding terms derived from 
the Schwarzschild geodesic equations (in the appropri- 
ate coordinate system). This enforced the Schwarzschild 
ISCO in the test-particle limit, but caused deviations 
from this value for finite-?7. Similarly, Ref. [28| intro- 
duced Pade approximants to improve the convergence of 
PN-based templates (in part by again enforcing agree- 
ment with the test-particle limit). The EOB formalism 
provides the most successful version of this idea by mod- 
eling the two-body dynamics in terms of a Hamiltonian 
that is based on a particle with reduced mass /1 = rjM 
moving in the "?7-dcformed" Schwarzschild background 
of a central mass M. It is in light of these resummation 
approaches that the ability of the Blanchet-Iyer ISCO 
condition to predict the Schwarzschild ISCO is surpris- 
ing (and perhaps not widely appreciated). 



A. Summary of results 

It is possible that the ability of the Blanchet-Iyer ISCO 
condition to predict the Schwarzschild ISCO is coinci- 
dental. One of the primary objectives of this study is to 



4 Indeed, one can see from Eq. I I1.2I I that the Schwarzschild ISCO 
frequency arises only from the 1PN equations of motion; the 2PN 
and 3PN terms affect only the 0(rj) corrections. Note also that 
in deriving this result, it was crucial to express Co in terms of 
the gauge-invariant observable x rather than a gauge-dependent 
radial coordinate |24H . 

5 The critical point defined in this way is sometimes called an 
ICO (innermost circular orbit). See Sec. II B of 21] (as well as 
Sec. IV A 2 of |25l |') for a discussion of the difference and rela- 
tionship between the ISCO and ICO. In the rest of this article, 
I will refer to both terms as an ISCO. 
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test this by extending the Blanchet-Iyer ISCO condition 
[Eq. (|1.2|) ] to the case of spinning (nonprecessing) bina- 
ries. This calculation is performed in Sec. HU The result 
is given by [see also Eq. (12.29)) below] 



Co = 1 - 6a; + x 3/2 f 14 



S £ 



6 



5m S' 



M V/ 2 ") 



IJ "- 



r 5/2 



(22 + 32? ? ) 



/397 



(5m E£ 
123 



1G 



(18 + 15?7) 
- Urf 



(1.5) 
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where S£ = £ • S c , Y,\ = £ ■ £ c , = 
unit vector along the direction of the Newtonian orbital 
angular momentum, S c = + S§, S c = M(S2/m 2 — 
Sf/mi), Sg = (1 + ma/mi)SS + (1 + mi/m 2 )S|, <5m = 
mi — m-2, and = Xa to a^5i are the individual spin 
angular momenta for body A — 1,2 with dimensionless 
spin parameters \a an d un it direction vectors s° A . This 
condition is derived from the 3PN equations of motion, 
including all explicitly known spin terms up to 2.5PN 
order. 

In the test-particle limit (77 — »• 0), the ISCO determined 
from Eq. (II .5|) can be compared with the ISCO of the 
Kerr spacetime 29]. This comparison can be performed 
by deriving a condition analogous to Co from the Kerr 
metric, expanding the result in powers of the BH spin 
Otf")) an d comparing to Eq. (|1.5p (see Sec.|Hl]for details). 
The resulting comparison shows that the two conditions 
agree up to the order to which the PN spin corrections 
are known. This comparison is also shown graphically 
in Fig. [TJ Note the large improvement in comparison 
with the 3PN energy function [which includes spin cor- 
rections; see Eq. (|Alj) ]. Presumably, if higher-order spin 
corrections in the PN equations of motion were included, 
the error in comparison with the Kerr ISCO for large 
values of \x%\ would improve. This excellent agreement 
suggests that the standard PN equations of motion are 
able to exactly recover some strong-field results. 

A second objective of this article is to calculate the 
shift in the ISCO frequency due to conservative effects 
(Sec. HVjl . In particular, two types of conservative effects 
are considered: the first due to the GSF, and the second 
due to the spin of the test-particle. As discussed above, 
the conservative GSF ISCO shift was computed in [3, HH 
for a Schwarzschild background, and compared with var- 
ious PN calculations in [21]. Here we focus on the ISCO 
shift in the Kerr background, for which GSF calculations 
are not currently available. Instead, we make predictions 
for what that ISCO shift might be according to two ana- 
lytic approaches: the ISCO condition in Eq. (|1.5[) above 
and the recently developed spinning-EOB formalism of 
Barausse and Buonanno [2(| (Sec. IIV Aj) . (In the latter 
case, the ISCO shift was calibrated to match the exact 
Schwarzschild result [II [l5l].]) 
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FIG. 1. (color online) . Comparison of three different methods 
for computing the ISCO of a nonspinning test-particle in the 
Kerr spacetime. The solid (black) curve (labeled Kerr) refers 
to the exact result for the Kerr ISCO [Ref. H<J or Eq. (|37Fj) 
here]. The dashed (red) curve (labeled Co) is the r\ — > limit 
of the gauge-invariant ISCO condition derived here [Eq. (|1.5[) 
or |[U5])]. The dash-dotted (blue) curve (labeled E3PN) is the 
ISCO computed by minimizing the 3PN circular-orbit energy 
[Eq. (|A1|) ]. The inset shows the fractional errors of the E3PN 
or Co curves with respect to the Kerr curve. 



To quantify these conservative ISCO shifts, we expand 
the ISCO frequency as [Eq. P~3)) below] 



GSF 



Mn = m 2 n K ( X 2)[i + vc 

+ 0( V 2 ) + O(xiv') 



(X2)+TOl^ USPm (X2) 



COspin^ 
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where fi K (x 2 ) is the Kerr ISCO frequency [Hj|. The shift 
in the ISCO due to the conservative GSF is parametrized 
by the function c GSF (%2)- In Sec. IIV 51 this function is 
calculated via the EOB and Co approaches; the results 
are presented graphically in Fig.[2]and tabulated in Table 
HJ It will be interesting to compare these numbers with 
future GSF calculations in Kerr. 

In Sec. IIV CI the function c COspin (x2) is also calcu- 
lated via the EOB and Co approaches^ However, in 
this case the EOB calculation via the Hamiltonian in [2(| 
yields the exact (fully relativistic) result. This is because 
this Hamiltonian reproduces the Papapetrou-Mathisson- 
Dixon equations of motion [33Tj4fj| in the small-?7 limit. 
An analysis of the ISCO shift directly using the Papa- 
petrou equations is also presented in Appendix [Bj the 
results are identical to those obtained from the EOB 
Hamiltonian (providing further confirmation of the work 



6 The ISCO for a spinning test-particle in Kerr was previously con- 
sidered in |3C| , but those authors focused on unphysically large 
values of the test-particle spin and did not explicitly compute 
the shift parameter c COs P in ( see also rjjrjjl]) 
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in USUI). These results are shown in Fig. [2] and Table Q] 
In the Schwarzschild case a fully analytic analysis of the 
Papapetrou equations is straightforward and presented 
in Appendix IB 31 It shows that the Boyer-Lindquist ra- 
dial coordinate of the Schwarzschild ISCO is shifted by 
O(ximi): 



6m 2 



m^Ximi + 0(xtmi), (1.7) 



and frequency shift of the ISCO due to the point-particle 
spin is given by 



COspin 



(0) = 



= 0.306186 



(1.8) 



Interestingly, this ISCO frequency shift is exactly repro- 
duced by the Cq ISCO condition, again showing that the 
standard PN equations of motion are able to exactly re- 
produce a strong-held result. (If \ 2 ^ 0, the exact result 
is only approximately reproduced by the Co condition 
because the PN spin terms are explicitly computed only 
to 2.5PN order; see Fig. [2) 

Section[V]discusses some conclusions of this study. Ap- 
pendix [5] compares the test-mass limits of several PN 
quantities (the orbital energy, angular momentum, and 
Keplerian relation) with the analogous quantities com- 
puted from the Kerr metric. 



A. Equations of motion and the relationship 
between spin variables 

The spin-orbit contributions to the equations of mo- 
tion given in [44j are expressed in terms of spin vectors 
SJ4 (A = 1,2) whose magnitudes y^ c m 2 A do not remain 
constant with time. (Note that Refs. [44l S3] do not use 
the superscripts "nc".) An alternative set of spin vari- 
ables are defined in Eq. (7.4) of [47| [also Eq. (2.21) 
of 46]] and have the property that their magnitudes are 
constant^ This choice of spin variables causes the spin- 
precession equations to take a convenient form and is 
generally preferred in computations. Here we denote the 
nonconstant- magnitude spin vectors of each body by S^ c , 
and the constant- magnitude spin vectors by S^. We also 
dehne the spin combinations 



S c = si 
E c = M 



^2 

m 2 



si 



(2.2a) 
(2.2b) 



and analogous relations for S nc and S nc . 

The relationship between (S nc ,£ nc ) and (S c ,£ c ) is 
given by Eqs. (2.22) of gH], 



2S nc + ^S nc 



M 



+ 0(c- i ), (2.3a) 



II. GAUGE-INVARIANT ISCO CONDITION 
FOR SPINNING BINARIES 



Followingthe stability analysis of the PN equations of 
motion in [24, 26j, we can generalize the gauge-invariant 
ISCO condition derived by Blanchet and Iyer [24j to the 
case of spinning, nonprecessing binaries. 

We begin by writing the conservative PN equations of 
motion for two spinning point-masses as 
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(2.3b) 



where r is the orbital separation in harmonic coordinates. 
The inverse relationship is given by 



gnc = gc 



1 



M 



-'I- 



2S C 



M 



v . S c + ^ v . S c 

M 



+ 0(c- 4 ), (2.4a) 



On the first line we list the nonspin terms to 3PN order 
(see [ll| for references); note that the radiation-reaction 
terms at 2.5PN and 3.5PN order are not present since 
we are only concerned with conservative corrections to 
the ISCO. The spin-orbit (SO) term at 1.5PN order and 
the spin-spin (SS) term at 2PN order were first derived 
in I42l| . The SO term at 2.5PN order was first derived in 
43]. Here I use the forms given in Eqs. (5.7) of [44]. The 
2PN order quadrupole-monopole (QM) term was derived 
in [45j]; Ref. [46[ shows how to concisely combine this 
term (when specialized to black holes) with the 2PN or- 
der spin-spin term [see their Eq. (3.8)]. 



5m 
~M 
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M 



S c + (1 - 2?/)£ c 



v ■ S c + (1 - 3?,)v • E c 



vj+0(c- 4 ), 



(2.4b) 



7 Throughout this section all of our spin variables are contravariant 
vectors. In [46ll these are denoted with an overbar. Note that the 
spin variables used in Kidder [4£l are the constant-magnitude, 
contravariant spin vectors denoted here. Note also that we 
use the notation S for the quantities denoted A in [46l . |4S(I . 
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where the powers of c were added to show that the correc- 
tions to the spins are a relative 1PN order effect. We also 
note the relationship between the individual spin vectors 




(2.5b) 

Since the spin variables differ at 1PN order, the equa- 
tions of motion (but not the equations of precession) will 
have the same form for the 1.5PN and 2PN spin terms 



(aside from the replacements S^ c -H- S^), but the 2.5PN 
and higher-order spin terms will differ depending on the 
choice of spin variables. Throughout this paper the su- 
perscripts "c" and "nc" are sometimes dropped where 
either index would be appropriate. 

The 2.5PN spin-orbit corrections to Eq. (|2.1[) are given 
in Eq. (5.7) of [44[ in terms of the variables S nc and 
S nc . The equivalent expressions in terms of the constant- 
magnitude spin variables are found by substitutin g th e 
relations (J2U) into the 1.5PN SO term [Eq. (5.7a) of [HI], 
and combining the result with the 2.5PN SO term in 
Eq. (5.7b) of HU (into which the substitutions S -» S c 
and S — > S c can be made since we only require accuracy 
to relative 2.5PN order in the spin terms). The resulting 
SO contributions to Eq. (|2.1[) in terms of the "c" spin 
variables are 



B 

so 



1.5PN 



A777 

12(S c ,77,t,) + 6 — (£ c ,n,«) 



f) 777 u 777 1 

9(nv)n x S c + 3 (nv)n x £ c - 7v x S c - 3 v xS c . (2.6a) 



B 

so 



2.5PN 



, ,0 9 M , A 5m, „ .( , , 9 , M ( 29 
(S c , n, v) [ -30-q(nv) 2 +2Ar]v 2 (44+25?/) J + — (£ c , 77, v) ( -15rj(nv f '+12rjv 2 ( 24 + — 77 



777 

(S c , n, w)(-9 + 977) + — (£ c , n, w)(-3 + 67,) 



5m ( M , „. 3 . 

— (4-77(77^) + -(777;) (t,S c ) 



+ (nv)n x S c 



+ -j7( ro ) n x S 



15?7(777;) 2 + I2777; 2 - 



M 







(-!•<); 


+ v x S c 



3 M 
-(nv)(vS c )(l - 77) - 8— 77(77^) 
2 r 



45 , x n n M , 

- -7TV(nv) + 21r]v 2 - 7— 4 + 3rf\ 
2 r 



33 M 

—T](nv) 2 H (24 + II77) - U777; 2 

2 r 



M 



-v x S c 



M 



9r)(nv) 2 - 7t]v 2 H (12-1 77) 



(2.6b) 



In the above equations we define additional notation fol- 
lowing [44]: the unit vector n = x/r points in the direc- 
tion of the relative separation vector x = yi — y2 ; v = x 
denotes the relative orbital velocity; scalar products of 
vectors are denoted by (ab) = a • b; and the mixed prod- 
uct of three vectors is denoted by (a, b, c) = a • (b X c). 

The sum of the spin-spin and quadrupole-monopolc 
terms is given in Eq. (3.8) of [46j|, 
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Note that Eq. (|2.7[) has the same form in terms of the 
"nc" spin variables; it is also only valid for Kerr BHs 
as the value for the Kerr quadrupole moment was used. 
Higher-order spin-spin corrections have recently been 



computed in Refs. [49-59], but the explicit equations of 
motion have not yet been derived. 



B. Restriction to the nonprecessing case 

Now we restrict to nonprecessing orbits in which the 
individual spin vectors are aligned or antialigned with 
the direction of the Newtonian orbital angular momen- 
tum vector I = Ln/|Ln|- We additionally define the unit 
vector A = £x n. Vectors can then be decomposed on the 
orthonormal basis {n, \,£} as in S = S n n + S\X + Si£; 
similar relations hold for S and So (in either spin repre- 
sentation), as well as for v. The restriction to nonprecess- 
ing orbits having a fixed orbital plane in the direction of 
£ then implies the following relations: 



•2 1 2-2 



v = + rtp\, 'i 
(nv) — r, n x v = rtp£, 

S = Sg£, X = ^i£, Sq = Sq (£, 



(2.9a) 
(2.9b) 
(2.9c) 
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(S,n,v) = rtpSi, (£, n, v) = rip^e, (2.9d) 

n x S = -SiX, n x S = -E e X, (2.9e) 

v x S = Si(rtpn — rX), v x S = T.e(r(pn — rA), 

(2.9f) 

(nS) = (nS) = (uS) = («E) = (n5 ) = 0. (2.9g) 



Eqs. (181)-(196) of Blanchet's review article gj. Denot- 
ing Blanchet's expressions by ,4b, ns and Sb,ns, ignoring 
the dissipative terms at 2.5PN and 3.5PN orders, and us- 
ing the form of the equations without the 3PN logarith- 
mic terms, the nonspin terms in Eqs. (|2.11|) are related 
to Blanchet's by 



The above relations allow the conservative PN two- 
body equations of motion to be put in the following form: 
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(2-10) The 1.5PN spin-orbit terms are found to be 
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The nonspin terms ,4. NS and B NS have been explicitly cal- 
culated by various authors. The results can be found in 



and have the same form in terms of the "nc" variables. 
The 2.5PN spin-orbit terms in both spin variables are 



A.5PN = — ( r( P) 



A 



SO.nc 
2.5PN 



M 

■—{rip) 
r 



— (20 + 1477) + f9 - y77 ) f 2 - 10t7(ty>) 2 
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one 
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M 2 
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M 
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Sm T, c e 

atm 2 

(2.14a) 
Sm Sf" 
AfM 2 

(2.14b) 
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M 
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„SO,nc _ M 
°2.5PN — ' 



M 
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r 

M 

-2— (2 + 317) + 7;f 2 + (9 - 2Ti)(ripY 
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CC 

M 2 

cnc 

M 2 



M 

—6 — 77 — ?7r 2 + (3 — rj)(ripY 
r 

M 

—4 — 77 — 77/ + (3 — rf)(rLpY 
r 



Finally, the spin-spin + quadrupole-monopole pieces are 

2 / CC \ 2 



, ss+qm _ 3 / M\ / 5g^\ 
^ 2PN ~ 2 I r ) \AP ) ' 



*2PN 

R SS+QM _ n 
°2PN — u ' 



(2.16a) 
(2.16b) 

and have the same form in terms of the "nc" variables. 



C. Perturbing the equations of motion 

Having simplified the equations of motion, we now wish 
to study perturbations about the circular orbit solutions. 
We first reexpress the equations explicitly in terms of the 
polar coordinates (r, if) of the relative position vector. 
Differentiating the expression for the velocity vector in 
Eq. (I2.9a[) and using h = ipX and A = — <pn, the compo- 



nents of Eq. (|2 . 10[) along n and A are given by 
V/ (l + ^ tot )+r^ 2 , 



1 f M 



r \ r 



$ = --[ ^B tot + 2f<p 



(2.17a) 
(2.17b) 



This system can be reexpressed in first-order form by 
defining u = r and uj = <p, resulting in three first-order 
equations in the variables (r, u,u>). 

Circular orbits correspond to the conditions r = u = 
Cj = 0. In particular, the condition u — and Eq. (|2.17a[) 
imply the following implicit relationship for the circular 
orbital frequency: 



Wo 2 = ^[l+X ot (^o,w )], 

r 



or, in terms of the PN parameter x = (Mwo) 2 / 3 , 
x = 7 [l + 4 ot (7^)] 1/3 , 



(2.18) 



(2.19) 



where a subscript refers to quantities evaluated along a 
circular orbit and we have defined another PN expansion 
parameter 7 = M/r. 

Equation (12.191) provides an implicit relationship be- 
tween the two PN expansion parameters 7 and x. Later, 
we shall need an explicit PN expansion for 7 in terms of x. 
To derive this relationship from (I2.19|) , we first substitute 
a 3PN series expansion with undetermined coefficients, 



7 = x(l + cix + ci. 5 x 3/2 + c 2 x z +c 2 .5X i,rz +c 3 x s ), (2.20) 



5/2 



into the right-hand-side of Eq. (|2.19p . Next we series 
expand the result in x to 3PN order, and equate the co- 
efficients of like powers of x on both sides of the equation. 
This results in a linear system of 5 equations for the 5 
unknowns in Eq. (|2.20l) . Solving this system easily yields 



1 + a; 1 



r 3/2 



,5/2 



+X° 



65 
8 

1 — 77 

3 9 ' 

2203 41 
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(2.21) 



In terms of the nonconstant spin-magnitude variables, 
the 2.5PN order term in the above equation should be 
replaced with [see Eq. (6.3) of [U 



r 5/2 



13 2 

y + 9 77 



one 

M 2 



8m Ej c 

Hm 2 



(2.22) 



while the 1.5PN and 2PN order spin terms have the same 
form with "c" replaced by "nc" . 

Now we examine linear perturbations to the equations 
of motion (I2.17[) about circular orbits parametrized by 
(ro,wo)- Introducing a small expansion parameter e, we 
substitute the following expansions into Eqs. (|2 . 1 7|) : 



r = r + eSr, 
u = + eSu, 
lu = uj Q + eduj, 

and linearize. In doing so we expand A tot 



(2.23a) 
(2.23b) 
(2.23c) 



MS 



A* 



\ tot 



dA b 



di 



Sr 



di 



Su 



OA" 



duj 



Sco, 



(2.24) 

and likewise for B tot . From the explicit form of A tot and 
B tot , one can verify that 



Q A tot 




Q B tot 




dB tot 


du 





dr 








0. 



(2.25) 



Then, at O(e ), the equations of motion reduce to 
Eq. f2TT81 and = 0. At C^e 1 ), we have the system 



Sr = Su, 



(2.26a) 



Su = a>oSr + (3qSu>, 
Slu = 7o^u, with 



a 



3wn - 



M dA t 



Or 



M dA tot 
Po = 2r w j— — 



1 / MdB tot 

70 = 2o; 

ro 



du 



(2.26b) 
(2.26c) 



(2.27a) 
(2.27b) 
(2.27c) 



where 70 is not related to the 7 = M/r defined earlier. 

Now we assume a perturbation of the form Sq = E q e M 
[where q = (r, u, 10)] and substitute into Eqs. (|2.26l) . 
resulting in a linear algebraic system for the E q and 
the eigenvalue A. A trivial solution corresponding to 
A = is E u = and E r = —(/3o/ao)E u ; this rep- 
resents a nonoscillatory displacement from one circu- 
lar orbit to another. The remaining eigenvalues are 
A = ±[— (ceo + /3o7o)] 1//2 - If the argument of the square- 
root is positive, then the resulting solutions are stable. 
The condition for the existence of stable circular orbits 
can therefore be expressed as 



C 



-ao - ft>7o > 0, 



(2.28) 



and the equality Cq = defines the ISCO. 

Using Eqs. ([2~27l) and (|2~TT|) . eliminating r via (j2~2"T]) . 

and expanding to the appropriate PN order, one can ex- 
press the stability condition explicitly in terms of x, yield- 
ing the following gauge-invariant condition for the ISCO: 



x 3 \ M 2 M M 2 



14?? -3 



CC 

°o.e 
W 2 



r 5/2 



/397 123 2 \ a " 



(2.29) 



In terms of the nonconstant-magnitude spin variables, 
the 2.5PN spin-orbit term is replaced with 



r 5/2 



one 

£ 

M 2 



(13 + 30r;) 



Srn 

M M 2 



(9 + 14??) 



(2.30) 



while the 1.5PN and 2PN spin terms have the same form 
with "c" relabeled to "nc" . Note that in the nonspinning 
case, Eq. (|2.29l) reduces to the 3PN gauge- invariant sta- 
bility condition of Blanchet and Iyer [24j [their Eq. (6.41) 
or Eq. (TOI) here] . 
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III. COMPARISON WITH THE KERR ISCO 

In the nonspinning case, Eq. (|2.29[) reduces in the test- 
mass limit to 



Cn = 1 - 6.T. 



(3.1) 



The PN ISCO criterion Co = in this case clearly repro- 
duces the exact Schwarzschild ISCO, x = 1/6. We wish 
to determine if Eq. (|2.29j) similarly reproduces the Kerr 
ISCO. 

Recall that the Kerr ISCO radius in Boyer-Lindquist 
coordinates is given by [29( 

3 + Z 2 - sign( X K)[( 3 _ Zl )(3 + Z 1 + 2Z 2 )\ 1 I' 1 , 



m 2 



Zx = 1 + [1 - (X^) 2 ] 1/3 [(l + Xf) 1/3 + (1 - X?) 1/3 ], 
Z 2 = [3( X f) 2 + Z 2 1 } 1 / 2 , (3.2) 

where the mass of the Kerr BH is denoted m2, and its 
dimensionless spin is X2 ^ [ — 1, 1] (with negative values 
corresponding to point-particles with retrograde orbital 
motion) An expression equivalent to Eq. (|3.2[) can be 
found by differentiating the reduced particle energy [29| 



E 



E 
mi 



1 - 2w Bh + xfw^ 2 
1 - 3ui B l + ZxfwBL 



(3.3) 



where wbl = m,2/rBL and tbl is the Boyer-Lindquist 
radial coordinate. Some simple algebraic manipulation 
of dE/dr B h = yields 

Cf = 1 - 6w BL + Sx^bl " Hxf) 2 w 2 Bh = 0. (3.4) 

Solving this equation for tbl produces results identical 
to Eq. (I3.2[) . But note that since u>bl depends on a co- 
ordinate radius, Eq. (|3.4p is clearly not a gauge-invariant 
expression. 

To derive a gauge- invariant version of Eq. (I3.4[) , wc first 
define the variable X = |rn 2 ^ K | 2 ^ 3 5 which is analogous 
to the PN parameter x (in the test- mass limit, x — >• X). 
The frequency fi K = dip/dt refers to the circular-orbit 
angular frequency seen by a distant observer and follows 
from the Kerr geodesic equations [Eq. (2.16) of |29j]: 



m 2 fJ K = sign(x 2 ) 



w 



3/2 
bl 



1 , K 3/2 ' 
1 + Xf W BL 



(3.5) 



Defining j3 = 1 — X2X 3 / 2 , we invert Eq. Q3.5P to obtain 
wbl = -jrpr- (3.6) 



In the notation of the previous section, the BH spin angular 
momentum is = x^ m 2^> where, in our restriction to non- 
precessing circular orbits, the orbital angular momentum points 
in the £ = z direction and we choose = z. 



Substituting this result into Eq. (|3.4[) . we arrive at the 
gauge-invariant relation 



r K = 1 - — 

^0 - 1 p 2 /3 



6 



xf 



X 1 ' 2 

01/3 



X 1 ' 2 
/3 1 / 3 



(3.7) 



For x% = we easily obtain the Schwarzschild value for 
the ISCO frequency (X = 1/6). One can verify numeri- 
cally that solving C^ = as a function of x^ reproduces 
the ISCO frequency computed from Eqs. (|3.2[) and (|3.5[) 
for all values of x^ £ [—1,1]. 

Now we wish to compare the test-mass limit of the 
ISCO condition derived in Eq. (|2.29[) with the gauge- 
invariant Kerr ISCO expression in Eq. (|3.7I) . Note that 
Eq. (|3.7|) is valid for arbitrary spin, while Eq. (I2.29|) is 
limited by the PN order to which spin terms have been 
computed in the equations of motion (currently 2.5PN 
order). To allow a meaningful comparison, we must ex- 
pand Eq. ()3.7|) in the spin parameter X2S yielding 



Cf = 1 - 6A + x^S^ 372 - 4A 5/2 ) 
+ (x^) 2 (-3A 2 + 8A 3 - 10A 4 /3) + 0[( X 



K\3l 



(3.8) 



We can also perform a PN expansion of Eq. Q3.7P in X, 
which results in 



,K^2 j^2 



C < = l-6X + 8 X fX 3 / 2 -3( X ^ 

- 4 X fx"/ 2 + 8( X f) 2 X 3 + 0[( X f) 3 Xy 2 } 



(3.9) 



Note that both expansions give consistent results at the 
appropriate orders in X2 an d X. This is especially inter- 
esting because in Eq. (|3.8[) . no PN expansion has been 
made. It also suggests the presence of additional self-spin 
terms at 3PN and 4PN orders in the equations of motion 
(in addition to the currently known 2PN-order terms). 
Equation (|3.9[) suggests that cubic self-spin interaction 
terms will not appear until 3.5PN order. 

The above expansions can now be compared with the 
test-mass limit of Eqs. (l2~^9"j) and (j2~3(Jl) . Taking 77 ->• 0, 
Sm/M -> -1, and {S e /M 2 ,^/M 2 , S 0/ /M 2 ) -> X 2, the 
result is 

C = l-6x + 8 X 2x 3/2 -3xW ~4 X 2X 5/2 + 0{x 3 ). (3.10) 

This is valid for either choice of spin variable [x 2 c or x 2 )- 
Comparing with Eq. (13.81) (and identifying X with x and 
xf" with X2), we see that the PN gauge- invariant ISCO 
condition (|3.10p agrees with the Kerr ISCO condition up 
to the PN order (2.5PN) to which we know the spin terms 
in the PN equations of motion. 

Figure [1] compares different methods for computing 
the ISCO frequency (in the test-mass limit): (i) the ex- 
act Kerr expression [computed from solving Eq. (|3.7p or 
plugging Eq. (|3.2p into Eq. (|3.5[) ]: (ii) solving the gauge- 
invariant ISCO condition in Eq. (l2li9"|) or (j3~10l) ; and (iii) 
finding the minimum of the PN circular-orbit energy with 
nonspin terms to 3PN order and spin terms to 2.5PN or- 
der [Eq. (|A1P ], [The ISCO frequencies for approaches 
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(i) and (ii) are also listed in Table HI] The method us- 
ing the gauge-invariant condition Co agrees exceptionally 
well for all spins up to X2 ^ 0.5. In the nonspinning case 
(X2 = 0), the agreement is exact. For nonzero spins, 
agreement with the exact Kerr result is limited by the 
fact that we only know the spin terms in the equations 
of motion to 2.5PN order. Note also that for small \x2\1 
the error is symmetric about \i — 0. This is in contrast 
with the ISCO computed from the 3PN energy function, 
for which the error increases (nearly) monotonically with 
increasing ISCO frequency (or decreasing radius). This 
indicates that the ISCO computed via C is limited not 
by finite-PN corrections but by finite-spin corrections. 

In Appendix [X] we examine how other PN expressions 
agree with their Kerr-spacetime counterparts. We find 
that test-mass limits of the circular-orbit energy and the 
Keplerian relation 7(2;) agree with their Kerr analogs if 
we identify Xsf with either choice of spin variable. How- 
ever, the PN orbital angular momentum only agrees with 
its Kerr analog if we identify X2" with x\- 

IV. CONSERVATIVE SHIFTS IN THE ISCO 

Consider the general behavior of the ISCO frequency 
when the test-particle has a non-negligible mass and spin 
(but assume that all spins are aligned or antialigned with 
the orbital angular momentum). The ISCO frequency 
can be split into the following pieces 

m 2 n = n K (x2)+SCl GSF ( X 2,q) 
+ SCl co ^( X 2,q, Xi) + ^ GSF+COspin (x2, Q, Xi), (4.1) 

where f2 K is the Kerr ISCO frequency in units of mi 
[given by Eqs. (J32J| and or Eq. ([XT])], 5Cl GSF and 

<5fi COspin are corrections to this frequency (also in units 
of m 2 ) due to the conservative GSF and the spin of the 
smaller compact object, and Jf2 GSF + COs P m jg a correction 
that results from cross-terms between both effects. If we 
assume that the mass ratio q = m,\jmi < 1 is small, then 
we can rewrite Eq. (|4.1[) as 

m 2 n = n K (x2)[i + qc' GSF (x2) + <zxic COspin (x 2 ) 

+ 0(q 2 ) + 0{ X iq 2 ) + 0{xW)]. (4.2) 
Multiplying by M/m 2 and using i] = q + 0(q 2 ) [l9| yields 

n = mq = n K ( X 2)[i + vc GSF (x2) + vxic COspin (x2) 

+ 0( V 2 ) + 0( X iV 2 ) + 0(xW)l (4-3) 

where c GSF = 1 + c' GSF was labeled c^ on in [H, EH for 
the X2 = case. 



In the remainder of this paper and unless stated otherwise, all of 
the spin variables refer to the constant-magnitude spins. 



In the remainder of this section, we shall concern our- 
selves with the calculation of the coefficients c G (X2) 
and c COspin (Y2) via the improved spinning-EOB Hamil- 
tonian of [2(| and the new gauge-invariant PN ISCO con- 
dition in Eqs. (|2.29l) . In particular, we note that the im- 
proved EOB Hamiltonian is constructed such that the 
coefficients c G Qg(0) and c , ^ 111 (X2) are exact. 

A. The improved effective-one-body Hamiltonian 
for spinning binaries 

Recently, Barausse and Buonanno have con- 

structed a new EOB Hamiltonian with the following fea- 
tures: (i) In the test-particle limit, the Hamiltonian re- 
duces to the exact Hamiltonian of a spinning test-body 
in the Kerr spacetime [40J (to linear order in the test- 
particle's spin; this limit of the EOB Hamiltonian pro- 
duces equations of motion and precession that are equiv- 
alent to the Papapetrou-Mathisson-Dixon equations (33r - 
[3§|). (ii) When PN-expanded, the EOB Hamiltonian re- 
produces the 2PN spin-spin and 1.5PN and 2.5PN spin- 
orbit couplings for arbitrary mass ratios, (iii) The Hamil- 
tonian includes an adjustable function K (rj) that appears 
in the spinning generalization of the effective metric func- 
tion A(r) [see Eqs. (6.9)-(6.11) of [2(|]; this function is 
adjusted to enforce agreement with the Barack-Sago con- 
servative GSF shift in the Schwarzschild ISCO 0, [H[ . 
[But note that this adjustment does not guarantee good 
agreement with the (yet uncalculated) conservative GSF 
shift in the Kerr ISCO.] (iv) For arbitrary mass ratios, 
this improved EOB Hamiltonian provides a well-defined 
prescription to compute the conservative two-body dy- 
namics and spin precession, (v) Finally, in the case of 
aligned or antialigned spins, this conservative dynamics 
produces a well-behaved ISCO for any mass ratio. 

The improved EOB Hamiltonian of [20( is compli- 
cated to write out explicitly. For the case of equa- 
torial (nonprecessing) orbits with spins aligned or an- 
tialigned with the orbital angular momentum, one can 
construct the Hamiltonian by starting with Eq. (6.1) of 
20] and carefully following their paper for the subse- 
quent chain of definitions (see Appendix C of [6(| for 
an alternate presentation). Once the EOB Hamilto- 
nian is constructed, the ISCO angular frequency can 
be computed from Eqs. (6.6)-(6.8) of [20j . Choosing 
units in which the total mass M = 1, I constructed 
a numerical code which computes the ISCO frequency 
^eob(?7, XI1X2) given the reduced mass-ratio 77, the spin 
of the test-particle xi> an d the BH spin X2- By con- 
struction, the resulting EOB ISCO has three important 
properties: (i) in the test-particle limit it reduces to the 

Kerr ISCO [fW(0, 0, X a) = ^ K (x 2 )]; (ii) in the non- 
spinning case it reproduces the exact conservative GSF 
ISCO shift [cggf^O) = c£ n fa 1.251]; and (iii) it correctly 
accounts for the conservative ISCO shift due to the test- 
particle's spin (this was explicitly verified in Appendix [Bl 
by directly analyzing the Papapetrou equations) . 
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TABLE I. ISCO quantities as a function of the dimensionless BH spin parameter \2- The second column denotes the standard 
Kerr ISCO angular frequency in units of m,2 [Eqs. (|3,2I) and (|3.5p ]. The third column is the test-particle limit of the ISCO 
frequency computed from the gauge-invariant ISCO condition Co [Eq. (|2.29[) or (|3.10[l ]. The fourth column is the conservative 
self-force ISCO shift parameter computed from the EOB ISCO frequency [Eq. (14. 4[) ]. The fifth column is the analogous quantity 
computed from the Co ISCO condition [Eq. (|4.5[) ]. The sixth column computes the ISCO shift parameter due to the spin of the 
test-particle (computed via the spinning-EOB ISCO frequency [Eq. I|4.6p ] , or directly from the Papapetrou equations [Appendix 
IB]). The seventh column is the analogous quantity computed via the Co ISCO condition [Eq. (|4.7[) ]. Note the perfect agreement 
of several of these quantities in the X2 — case, and the closeness in their values for small \2 < 0.6 (see also Figs. [1] and [2]). 



X2 


^isco 




GSF 
c EOB 


GSF 


COspin 
C EOB 


COspin 

c PN 


-0.99 


0.038 635 


0.038 015 


0.9486 


1.1903 


0.2313 


0.1945 


-0.9 


0.040 261 


0.039 681 


0.9449 


1.1961 


0.2364 


0.2020 


-0.8 


0.042 223 


0.041 694 


0.9423 


1.2043 


0.2424 


0.2110 


-0.7 


0.044 372 


0.043 901 


0.9422 


1.2148 


0.2487 


0.2205 


-0.6 


0.046 736 


0.046 331 


0.9458 


1.2282 


0.2553 


0.2308 


-0.5 


0.049 348 


0.049 016 


0.9550 


1.2453 


0.2625 


0.2417 


-0.4 


0.052 251 


0.051998 


0.9726 


1.2670 


0.2700 


0.2534 


-0.3 


0.055 496 


0.055 325 


1.0027 


1.2948 


0.2782 


0.2657 


-0.2 


0.059 149 


0.059 057 


1.0517 


1.3303 


0.2868 


0.2788 


-0.1 


0.063 295 


0.063 266 


1.1295 


1.3759 


0.2962 


0.2923 


0.0 


0.068 041 


0.068 041 


1.2513 


1.4349 


0.3062 


0.3062 


0.1 


0.073 536 


0.073 492 


1.4418 


1.5116 


0.3170 


0.3199 


0.2 


0.079 979 


0.079 750 


1.7400 


1.6118 


0.3287 


0.3328 


0.3 


0.087652 


0.086 978 


2.2072 


1.7434 


0.3414 


0.3435 


0.4 


0.096 974 


0.095 365 


2.9338 


1.9167 


0.3551 


0.3502 


0.5 


0.108 588 


0.105 125 


4.0204 


2.1441 


0.3699 


0.3499 


0.6 


0.123 568 


0.116470 


5.4310 


2.4388 


0.3856 


0.3382 


0.7 


0.143 879 


0.129 564 


6.3967 


2.8098 


0.4014 


0.3097 


0.8 


0.173 747 


0.144 421 


4.2785 


3.2524 


0.4150 


0.2590 


0.9 


0.225 442 


0.160 767 


-3.3671 


3.7337 


0.4152 


0.1837 


0.99 


0.364 410 


0.176 197 


-23.763 


4.1440 


0.2937 


0.0983 



EOB and PN predictions for the conservative 
self-force ISCO shift in Kerr 



The conservative self-force ISCO shift parameter de- 
noted c GSF in Eq. (14.31) is an especially interesting quan- 
tity because it is a gauge-invariant that can be calculated 
from self-force calculations. Barack and Sago [H], [H| 
have computed this quantity in the case of Schwarzschild, 
and in Ref. [21] this result was compared with multiple 
PN-based computations of the ISCO shift0 Gravita- 
tional self-force results are not yet available for the Kerr 
spacetime, but here we explore the predictions for the 
conservative GSF ISCO shift in Kerr given by two PN- 
based calculations: the spinning-EOB approach [2(| and 
the ISCO computed via the gauge-invariant PN ISCO 
condition Co [Eq. (|2 . 29[) ] . Based on the comparison study 
in [2l| , these two methods are the most viable approaches 
for computing the ISCO in the small-mass-ratio limit. 

Using the EOB ISCO frequency calculated from [20| 
as described above, the corresponding conservative GSF 



ISCO shift parameter can be computed via 



For other comparisons of PN and GSF results, see fsL Il9t [6ll463(l . 



c£ob(X2) = lim 



1 



r/->0 r\ 



^EOBfaO,X2) 
^ K (X2) 



(4.4) 



In the PN case a function ClcoiV) Xi> Xa) is computed by 
solving for the root of Eq. (|2.29j) numerically. The result- 
ing conservative GSF ISCO shift parameter is defined by 



„GSF 

-Co 



(Xa) = I™ - 

r;-s-0 T] 



n Co (o,o,x 2 ) 



(4.5) 



Note that in this equation the denominator contains the 
function Clc (0, 0, X2) rather than (l K (x2)- This is be- 
cause the gauge-invariant PN ISCO (lc does not reduce 
precisely to the Kerr ISCO (although it is very close for 
small to moderate values of \2\ see Fig. [1] and Sec. IIIip . 

The resulting values for CgQg(x2) and c < g SF (x2) are 
listed in Tableland plotted in the left- half of Fig.[l Note 
that while the EOB curve is calibrated to the exact result 
in the nonspinning case, there is no expectation that it 
will also predict the correct ISCO shift in the spinning- 
case. The function K(rj) will presumably need to be re- 
calibrated when GSF results for the Kerr ISCO shift are 
available. To further explore the behavior of Ceqb(X2)j 
have varied the value of K from to 4. Figure [5] shows 
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FIG. 2. (color online). ISCO shift parameters computed via the improved spinning-EOB Hamiltonian of [20| and the gauge- 
invariant ISCO condition in Eq. (|2.29[) . The left plot shows the ISCO shift due to the conservative gravitational self-force 
(GSF) as a function of the big BH spin \2 (the test-particle is assumed to be nonspinning in this case). The solid (blue) "EOB" 
curve uses the Hamiltonian from [20| [which is fit to the exact Barack-Sago (BS) result in the nonspinning case] and Eq. (14.41) . 
The dotted (green) "EOB (uncalibrated)" curve also uses this Hamiltonian, but the adjustable function is set to K{rf) = 1/2. 
The dashed (red) curve labeled "Co" is from Eqs. (|2.29[) and (|4.5|) . The right plot shows the ISCO shift due to the spin of the 
orbiting test-mass. In this case the "EOB" curve [Eq. (|4.6I) ] exactly reproduces the ISCO shift computed from the Papapetrou 
equations (see, e.g., Appendix iBl the Hamiltonian in [20( was constructed with this property). The "Co" curve [Eq. (14.711 ] 
agrees precisely with the exact result in the \2 = case. The difference between the c COspm curves for nonzero \2 arises from 
our limited knowledge of higher-order PN spin corrections. 



one of these "uncalibrated" choices \K{rj) = 1/2]. Vary- 
ing K over this range changes the location of the "peak" 
°f c eob(X2)- While the Barack-Sago result is no longer 
reproduced for other choices of K (the difference with 
the Barack-Sago value at X2 = gets especially large 
for K > 2), it is interesting to note that both the cali- 
brated and uncalibrated curves approach similar values 
when X2 — > ±1- 

It will be very interesting to compare future GSF cal- 
culations of the ISCO shift in Kerr with the results shown 
here. Strictly speaking, the values for c§^ F (x2) cannot 
be precisely compared with the "exact" \2 ^ GSF re- 
sults because the ISCO frequency in this case does not 
reduce precisely to the Kerr value. Still, for a large range 
of \2 (as quantified in Fig. [T]), an accurate comparison 
with future exact GSF results should still be possible. 
Note, in particular, that all three curves in the left-half 
of Fig. [5] roughly agree for \2 ^$ 0.2. This is perhaps 
indicative that the exact GSF results will lie near those 
values. These predictions are likely to be most accurate 
for \2 ~ — 1; varying K from to 4 near this value indi- 
cates cfg|(-l) 0.8-1.1. 



C. Conservative ISCO shift due to the 
test-particle's spin 



the point-particle. Using the EOB ISCO frequency, this 
quantity is calculated via 



COspin/ \ v 

"bob (X2) = hm 

r\— >0 



^EOBp7,Xl,X2) - »eob(t7,0,X2) 
VXi& k (X2) 



(4.6) 

Although the quantity c^q^ above is not exact (ex- 
cept for xi = 0), in this case the EOB Hamiltonian is 
constructed such that c EO g Pm (x2) is in fact the "true" 
value that would result from a calculation based on the 
Papapetrou-Mathisson-Dixon [33M40l j equations of mo- 
tionFH This was verified by an explicit calculation di- 
rectly based on the Papapetrou equations (Appendix |B|); 
the two methods give identical results for c pin (x2)- 
In the case of the Cq ISCO condition, 
compact-object spin ISCO shift via 



we define the 



COspin 
"Co 



(X2 



lim 

?7->0 



^Co (>7, XI1X2)- ^c (jh 0, X2) 
?7Xi^Co(0,0,X2) 



(4.7) 

where again the expression differs from Eq. (|4.6p because 
Clc (0, 0, X2 7^ 0) does not reduce to the exact Kerr ISCO. 

\X2) and c£° spm (x2) 



The resulting values for c^qb 111 ' 



are listed in Table Q] and plotted in the right-half of Fig. [2] 



It is also interesting to examine the ISCO shift pa- 
rameter c COs P in [Eq. (|4. 3p ] originating from the spin of 



11 Note that CgQg P '° does not depend on c^qj or the choice of the 
adjustable function K(r/). 
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Note, in particular, that in the Schwarzschild case the 
values for CgQg pm and c^? spin agree precisely with each 
other and with the analytic calculation in Appendix IB 31 

CE§B in (0) = cc° spin («) = ^Y = °- 306 186 ... . (4.8) 

This is a remarkable result. It indicates that the gauge- 
invariant ISCO condition Co not only predicts (i) the ex- 
act test-particle ISCO in the Schwarzschild case [24| , and 
(ii) the spin-expansion of the exact Kerr ISCO (Sec. IIIip . 
but it also predicts the exact shift in the Schwarzschild 
ISCO caused by the test-particle's spin. This shift is em- 
bodied in the (fully relativistic) Papapetrou-Mathisson- 
Dixon equations of motion, and it is rather unexpected 
that this shift could be predicted from an analysis based 
on the standard (nonresummed) PN equations of motion. 
Along with the other qualities mentioned above [and the 
closeness of c^ F (0) to the exact Barack-Sago result], this 
further indicates that there is a special quality to the 
gauge-invariant ISCO condition in Eq. (|2.29[) . 

In the spinning case, we see from Fig.[3]that c^° spm (0) 
starts to deviate from the exact result as \x%\ increases. 
This is due to the fact that the gauge- invariant ISCO 
condition Co is limited by the number of known spin 
corrections in the PN equations of motion. Once higher- 
order spin effects have been calculated and incorporated 
into these calculations, it is expected that the curves la- 
beled "C " in Fig. [T]and the right-half of Fig. [2] will even 
more closely approximate the exact results. 

V. CONCLUSIONS 



One of the most significant results of this study is 
that the PN ISCO condition in Eq. (gH5J) — in addition 
to reproducing the Kerr ISCO for small spin and the 
Schwarzschild conservative GSF ISCO shift with good 
accuracy — also exactly reproduces the ISCO shift due to 
the spin of the test-mass. (This agreement is truly exact 
only in Schwarzschild since the spin corrections in Co are 
only known to quadratic order.) This provides further 
evidence that the ability of the Co ISCO condition to 
predict strong-field results is not coincidental. However, 
it is somewhat mysterious as to why this ISCO condition 
is able to accurately predict these strong-field effects. 

In addition to explaining this agreement, future work 
could involve extending this study to more general orbits 
(such as precessing or eccentric binaries). The resulting 
conditions for the last stable orbit could then be com- 
pared with exact results from the Kerr spacetime. 

ACKNOWLEDGMENTS 

I gratefully acknowledge Luc Blanchet and Alessandra 
Buonanno for their helpful comments on this manuscript. 
For several useful discussions I also thank the members of 
the relativity groups at JPL and Caltech, as well as the 
participants of the "Theory meets data analysis at com- 
parable and extreme mass ratios" conference (Perime- 
ter Institute, June 2010). This research was supported 
through an appointment to the NASA Postdoctoral Pro- 
gram at the Jet Propulsion Laboratory, administered by 
Oak Ridge Associated Universities through a contract 
with NASA. 



The primary objective of this study was the exten- 
sion of the Blanchet-Iyer |24j] ISCO condition to the case 
of spinning, nonprecessing binaries [Eq. (|2.29[) ]. When 
the test-mass limit of this condition is compared with 
the exact Kerr ISCO, they are found to agree up to the 
order to which the PN spin terms are explicitly known 
[cf. Eqs. and (pTTUl) . and see Fig. [J. In addition, 

the conservative ISCO shifts were also computed using 
this ISCO condition and the spinning-EOB Hamiltonian 
of [H [see Table Hand Fig. [2]. 

The ISCO shift due to the conservative gravitational 
self-force should eventually be compared with the exact 
results from self-force calculations. This will allow an 
extension of the study in [2Ljj to the Kerr case, and will 
provide insight into the relative accuracies of the EOB 
formalism and the standard PN equations of motion. For 
example, in [2Ljj it was found that the Blanchet-Iyer ISCO 
condition more accurately reproduces the Barack-Sago 
ISCO shift than uncalibrated EOB methods. This excel- 
lent agreement in the Schwarzschild case could be coin- 
cidental, but it would be hard to dismiss if it were also 
true in the Kerr case. Comparison with exact self-force 
results in Kerr would clarify if the standard PN equations 
of motion or the (uncalibrated) EOB approach can more 
accurately predict strong-field, finite-?] effects. 



Appendix A: COMPARING THE PN AND KERR 
EXPRESSIONS FOR THE ENERGY, ANGULAR 
MOMENTUM, AND KEPLER RELATION 

In this appendix we examine the test-mass limit of var- 
ious PN expressions, and compare them with the equiv- 
alent expressions derived from the Kerr metric. In par- 
ticular, we wish to check if PN expressions using two dif- 
ferent choices for the spin variables reduce to the same 
Kerr result in the test-mass limit. 

The energy for circular, nonprecessing orbits is 




where the first two lines contain the nonspin terms [64 , 
|65| . the third line contains the 1.5PN spin-orbit term 
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[44 1471 . |48| and the combine d sp in-spin + quadrupole 
monopole term (for BHs only) [45l |46|. |48| . and the fourth 
line contains the 2.5PN spin-orbit term [44ll47j. In terms 
of the nonconstant-magnitude spin variables, the 2.5PN 
spin-orbit term can be written as [47| 



x 5 / 2 
M 2 



Sm„ r 



, (A2) 



while the 1.5PN and 2PN spin terms keep the same form 
but with "c" replaced by "nc". In the test-mass limit 
E FN reduces to 



£ FN (f}) 1 



mi 



;X2X 



3/2 



27 o 



r/ 9 675 

64 



(A3) 



where \i can be either x 2 ° or X 2 - 

The total energy of a point-mass in the Kerr spacetimc 
is given in terms of u>bl in Eq. (13.3[) . [Recall that E 
includes the particle's rest mass, so the orbital energy is 
E — 1.] Substituting Eq. (I3.6[) and expanding in X yields 



E-l 



X 



-x + -xfx 3 / 2 - X 2 

4 3 A2 



27 



(X 



K\2 
2 I 



■ s x fx 5/2 + X s 



675 _ 65 
~M ~ 18 



ix^f 



+27 x fX J / 2 + 0(X 4 )}, (A4) 



which agrees with Eq. (|A3j) to the expected order. 

The orbital angular momentum (specialized to equa- 
torial orbits) is given by Eqs. (6.10) and (7.10) of [47j], 



L £ = 



rjM 2 
-ZT/2 



l + x[- + - 
\2 6 

3/2 

■3 i 



X 



M 2 



35 
~6~ 



a , 27 197? t, 2 
x \ 

8 24 



5 8m, 



2 M 



+ 



M 2 



77 427 \ / 21 35 \ Sm 



(A5) 



Also note that the spin-orbit terms in L differ even at 
1.5PN order when one switches spin variable^ In the 
test-mass limit, these expressions reduce to 



hi 



m 2 



m\ x 

„2 



1/2 



3 10 

2 X Y X2X 



3/2 



L ■ £ 



27 n c 

Y + QM 



(Xlf 



7 X c 2 x^ 2 + 0(x 3 ) , (A7) 



m 2 



mi x 1 / 2 
27 



-X 



-X 



2 



;X?* 3/2 



- yX2 C * 5/2 + 0(* 3 )}, (A8) 



where the constants Cq M and Cq^ have not been explic- 
itly computed. 

The orbital angular momentum of a test-particle in 
Kerr is [29| 



L 

mi 



sign( X f)m 2 [1 - 2 X fw 3 B ^ + (xf?w 2 



BLJ 



1 - 3w B L + 2xf w BL 



Substituting Eq. (|3.6|) and expanding in X yields 



(A9) 



L 



X 1 



sign(x 2 < ) m 2 f- 3,, 10 



X 1 / 2 

? + (x^) 2 



2 X -^ 2XV2 



7 X fX 5/2 + X 3 



135 26. j 



K\2 



-^X^ 7/2 + 0(X 4 )}. (A10) 

Here we see that the Kerr angular momentum agrees with 
the test-mass limit of the PN expression only if we iden- 
tify X% w ith xl- Note also that Eq. (|A10j) provides the 
test-mass limit of the previously unknown 2PN and 3PN 
pieces of Eq. (|A5[) . 

We also check for agreement between the PN and Kerr 
versions of the Keplerian relationship (see also Appendix 
B of 01). The PN relation is given in Eqs. (|23Tj) and 
(I2.22[) . In the test-mass limit it reduces to 



with the last two lines replaced by the following expres- 
sion in terms of the "nc" spin variables: 



+ 



^3/2 
M 2 



23 3(5to 
~6 be "2M L< 



M 2 



77 259 
"~8~ + ~72 



nc 



33 Sm ync 



(A6) 



Note that the 2PN spin(l)-spin(2) term is zero HU, but 
the 2PN quadrupole-monopole contribution has not been 
computed. The 3PN nonspin terms are given in general 
form in [66j . but have not been specified to circular orbits. 



7712 



;X2X 



3/2 



x 2 [ 1 



9 

X-2 



:X2X 



5 / 2 +x 3 + 0{ X 2 2 x 3 ) 



(All) 



12 The total angular momentum J = L + Si/c + S2/cisa constant 
vector (up to 2PN order) that does not depend on the choice 
of spin variable. Since the individual spins contribute a 0.5PN 
correction to the total angular momentum, the 1PN corrections 
in the relations between spin variables [Eqs. 12.5H shift some 
terms into (or out of) the 1.5PN piece of L. 
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where xi can be either spin variable. Note that the PN 
radial coordinate used in the main text refers to harmonic 
coordinates (here denoted r#). To derive the Kerr-analog 
of this expression we first need the relationship between 
Boyer-Lindquist and harmonic coordinates |43l . l67j ] , 



x H 



iyn = (rBL - m,2 + ixf m 2)e lip sin6> B L, (A12a) 
zh = (rBL - ma) cos# B l, (A12b) 



H 



■A 



z h = ( r BL - rn 2 ) 2 



+ (X 2 V 2 ) 2 sin 2 : 



BL- 



(A12c) 



Specializing to the equatorial plane (9bl = 7r /2) and 
defining wh = rn 2 /rH, we have the relationship 



w H 



wbl ) 



(x>bl) 2 



(A13) 



Substituting Eq. (|3.6I) for wbl and series expanding in X 
yields 



w H = X \ 1 + X + -x 2 X 3 ' 2 + X 2 



(xf) 2 



1 - 



\xW 2 



1 - — (y k ) 2 
18 [X2 ' 



+xy 2 [2 x f-( x ff]+0(X±)}, (A14) 
which agrees with Eq. (|A1 1|) to the expected order. 



Appendix B: SPINNING TEST-PARTICLE ISCO 
SHIFT DERIVED FROM THE PAPAPETROU- 
MATHISSON-DIXON EQUATIONS 

In this appendix I discuss how to compute the ISCO 
for a spinning test-particle directly from the Papapetrou- 
Mathisson-Dixon equations (rather than from the EOB 
formalism of [20[). The results for the ISCO shift pa- 
rameter derived below agree exactly with the results dis- 
cussed in Sec. HVCl Table H and the right plot of Fig. [2] 
This provides further confirmation of the validity of the 
Hamiltonian derived in pol . 140} . Previous examinations 
of the ISCO of a spinning test-particle are given in [30ll3l1 | 
(see also [Hj]). The results here are more explicit, exact 
numerical values are given (Table IJ), and a fully analytic 
examination in Schwarzschild is presented. 



1. Papapetrou-Mathisson-Dixon equations for 
equatorial, nonprecessing orbits 



a form similar to the Kerr geodesic equations [Eqs. (2.19)- 
(2.25) of [HI]: 



dt 



S S A S -^- =a I 1 , 



3m2S z 



dy 

'd^ 



3m2S 2 
rE. 



J z - (a + s)E 



J z - (a + s)E 



A 



A' 



dv i — 
E S A S — = ±s/R s , where 
ar 



r 2 f i _ H^l 



A, = 1 



3m 2 s 2 r[J z - (a + s)Ef 

e! 



(Bla) 

(Bib) 
(Blc) 

(Bid) 
(Ble) 
(Blf) 



(r 2 +a 2 ) + a S (l + ^)]E-(a + S ^)j z , 

(Big) 

A = r 2 - 2m 2 r + a 2 , (Blh) 

where (i, r, 9, ip) are Boyer-Lindquist coordinated, r is 
the particle's proper time, and the conserved energy 
E = E/mi and total angular momentum J z = J z jm\ 
are given in Eqs. (2.10) of [68} . 

Note that the function R s can be rewritten in the form 

R s = B(r)[E - Ei(r, J z )][E — E 2 {r, J z )], (B2) 

where the roots Eyoof R s — are found by solving [see 
also Eq. (2.26) of [68j]] 



R s = P 2 - A {-§■ + [J z - (a + s)Ef 



P.* = 



aE 2 - 2j3E + 7 = 0, with 



m 2 



(r 2 + a z ) + as 1 + — - A (a + s) 



^{( a + S ^)[(, 2 + fl 2 ) +as (l + ^ 



(B3) 



(B4) 



(B5) 



7 = « + Jl-^ 



1 



(B6) 

Here (a, /?, 7) are not to be confused with any quantities 
defined earlier in this paper. The solution 



Ex = V c 



off 



/3 + v//3 2 - a 7 



(B7) 



Saijo et al. (6g| have explicitly derived the equations 
of motion of a spinning particle in the equatorial plane 
(6 = tt/2) of a Kerr BH. For a particle with spin angular 
momentum Si = sm\Z aligned with the BH's spin (S 2 = 
am 2 z) and the orbital angular momentum L z , the spin 
vectors remain constant and the equations of motion take 



13 Note that in the rest of this paper r denotes the harmonic radial 
coordinate. Also, to maintain some notational consistency, I con- 
tinue to denote the central BH mass by 1712 and the test-particle's 
mass by m\ ; in most of the literature on the Papapetrou equa- 
tions these quantities are denoted M and fi respectively. 
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corresponds to an effective potential for the particle mo- 
tiono Here we have taken the positive square root to 
ensure that the particle energy E = E\ — > 1 when r — > oo 
(in contrast to the negative root, for which E2 — > —1). 
This allows us to rewrite the equation for the radial mo- 
tion in the form 



r 2 - A(r,E, J z )[E-V cS (r, J z )} 



(B8) 



where, for this appendix only, an overdot means d/dr. 
The explicit forms for A and B can be inferred from the 
above equations but are not needed for the remainder of 
the analysis. 



General solution for the ISCO of a spinning 
particle 



as initial guesses for the solution. The resulting val- 
ues (7*0, Jo) are then used to determine the ISCO energy 
Eq = V e s(ro, Jo). The ISCO angular frequency is then 
found by substituting these quantities into 



n 



dip/dr 
dt/dr 



(B13) 



using Eqs. (|Bla[) and (|Blb[) . This procedure allows 
the ISCO frequency to be computed as a function of 

The ISCO shift parameter c spm (x 2 ) is computed as 
in Eq. (|4.6p . Note that in this case c GSF evaluates to 
zero (as expected) and converting variables from (q, m 2 f2) 
to (r),MQ) does not affect the value of c COspin (x 2 ) [see 
Eq. (|4.3p]. The resulting values for c COspin (x 2 ) are iden- 
tical to those listed in Table U under CEOB Pm (X2)- 



The conditions for circular orbits (defined as orbits 
with constant r) are that both r and r vanish. By differ- 
entiating Eq. (|B8[) and dividing by f, 



1 

r = — 
2 



- OA dV cS 

(E-V eS )—-A— 



(B9) 



we see that the conditions for circular orbits are equiva- 
lent to 



E = V e s(r,J z ) and 



dV c g(r,J z ) 
dr 



= 0. 



(BIO) 



To ensure that circular orbits are stable, we require that 
under a small radial perturbation of a circular orbit, 
7*0 + dr, the particle is accelerated back to its ini- 

tial configuration. Such a condition is equivalent to de- 
manding that the perturbed coordinate acceleration sat- 



isfy Sr 



-LinSr with cjg > 0, where cDq is the radial os- 



cillation frequency about the unperturbed orbit tq [this 
is equivalent to the analysis in Eqs. (|2.26p - (|2.28l) above]. 
In this case uiq is found by linearizing Eq. (|B9|) about the 
circular orbit tq. Computing dr/dr and evaluating along 
the unperturbed circular orbit yields 



Or 
dr 



A d 2 V cS 
2 dr 2 ' 



(Bll) 



The ISCO is found from the equality ojq = (note that 
A is nonzero for physically relevant parameter values). 

To evaluate the ISCO frequency, we first solve the al- 
gebraic system of equations 



dV cS (rJz) _ Q 
dr 



and 



d 2 V cS (r, J z ) _ 



dr 2 



= 



(B12) 



for the ISCO values of (r, J z ). This is done numerically, 
specifying a — X2 m 2, s — xiQ m 2, ^2 = 1, and us- 
ing r = rf sco [Eq. and J z = L(r* co ) [Eq. W 



3. Analytic analysis of the ISCO in the 
Schwarzschild, small-spin limit 

It is instructive to reexamine the above analysis of the 
ISCO, specializing to Schwarzschild (a = 0) and small 
spin (s/7712 *C 1)Ej Keeping terms linear in s, Eqs. (|B1[) 
reduce to 



dt 



E 



sm 2 J z 



Y _ 2ro 2 r 3 H _ 2m 2 \ 
r V r ) 

J Z sE 



0(s 2 ), 



E 2 



+ 0(s 2 ), 



-lV c t w (rJ z )} 2 - 



where (V^f 



2s 



EJ Z 



3m 2 




(B14a) 

(B14b) 

+ 0(s 2 ), 
(B14c) 

(B15) 



is the effective potential for Schwarzschild. Setting r 2 = 
yields a quadratic equation for E which, when solved and 
expanded in s, yields 



E = V 



schw, spin 



cir 



^schw 



sJ z 

r 2 



1 



3ms 



0{s 2 



r* \ r 

(B16) 

Note that this is equivalent to the O(s) expansion of 
Eq. (JB7) with a = 0. 

For circular orbits we solve the condition 
for J z = Jf hw + s8J z + 0(s 2 ), 



ot /-schw, spin / ^ 

dV eS /dr 
yielding the angular momentum for circular orbits, 

rJrri2 



^ycirc 



y/r — 3?TT-2 



s (r-2m 2 )(2r-9m 2 ) 2 

2 V?(r - 3m 2 )3/ 2 {S h 

(B17) 



Eq. (2.27) of [6S( has the wrong sign in front of the 07 term. 



15 Since s/m,2 = Xi<?i the small-spin limit is still quite accurate for 
EMRIs since g<l even if xi ~ !• 
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Substituting into Eq. (|B16|) and expanding in s yields the 
energy along circular orbits, 



^circ 



m 2 



r — 2rri2 



yjr(r - 3ra 2 ) 2r\r-2mi 



3/2 



+ o( s 2 ). 



(B18) 

To determine the ISCO we compute 
d 2 V^ hw ' spin /dr 2 = 0, substitute Eq. flBTZ]) for J 2 , 
expand to O(s), and solve for r = 6m 2 + sSf + 0(s 2 ). 
The resulting ISCO radius is 



M, m = 6m 2 — 2s\j - + 0(s 2 ). 



(B19) 



Substituting this result into Eqs. (|B18|) and (|B17[) gives 
the energy and angular momentum at the ISCO, 



To compute the ISCO frequency we expand Eq. (|B13I) , 

= j± / _ 2m,\ _ s± -2m*) ( rnjj\ 2 
r 2 E\ r J r 3 \ r 3£2 J V " 

(B21) 



substitute Eqs. (|B19|t and (|B20|) . and expand to O(s), 

1 s 



m 2 n isco = 6~ 3/2 + — — + 0{s 2 ). 

48 m 2 



(B22) 



Multiplying by M/m 2 and using s = XiQ m 2 and g = 
77 + 0(r] 2 ), we can write the shift in the ISCO as 

Mn iBCO = 6- 3 / 2 [l + r, + c™T n XiV + 0( V 2 )}, (B23) 

where the 0{rf) term would combine with the conserva- 
tive GSF shift (not computed here), and 



COspin 

C 



= 0.306186 217 847. 



(B24) 



Jl sco = 2V3m 2 + + 0{s 2 ). (B20b) 



This agrees with the Schwarzschild value found from 
three separate calculations via the EOB Hamilto- 
nian [Eq. (|4.6|) J. the gauge-invariant ISCO condition 
[Eq. (|4.7I) ]. and the numerical evaluation of the Papa- 
petrou equations (Appendix IB 21 above) . 
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